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Abstract. We present a GRANAT/SIGMA observation of the 
soft gamma-ray afterglow immediately after GRB 920723. The 
main burst is very bright. After ~ 6 s, the burst light curve 
makes a smooth transition into an afterglow where flux decays 
as f""'^. The power-law decay lasts for at least 1000s; beyond 
this time, the afterglow emission is lost in the background fluc- 
tuations. At least ~20% of main burst energy is emitted in the 
afterglow. At approximately ^^ 6 s after the trigger, we also ob- 
serve an abrupt change in the burst spectrum. At f < 6s, the 
ratio of 8-20 and 75-200 keV fluxes corresponds to the power 
law spectral index a = 0.0-0.3. At f = 6s, the value of a in- 
creases to a « 1 and stays at this level afterwards. The observed 
afterglow characteristics are discussed in connection with the 
relativistic fireball model of gamma-ray bursts. 
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1. Introduction 

Fast and accurate localizations of gamma-ray bursts by Beppo- 
SAX helped to establish the connection of GRB with the 
sources of decaying X-ray, optical, and radio emission (e.g. 
Costa et al. 1997, Van Paradijs et al. 1997, Frail et al. 1997). 
X-ray afterglows were found in 15 of 19 well-localized bursts; 
in most cases, the X-ray flux decayed as a power law of time, 
r^, with /3 ranging from -1.1 (GRB 970508, Piro et al. 1998) 
to -1.57 (GRB 970402, Nicastro et al. 1998). The power law 
decay of flux also is observed in the optical (e.g. Wijers et al. 
1997, Sokolov et al. 1998). This is a characteristic prediction of 
the relativistic fireball model of GRB (Meszaros, Rees 1993, 
Meszaros 1997, Waxman 1997, Sari et al. 1998). Indeed, the 
energy release in some GRB is enormous and sufficient to 
power the relativistic fireball (KuUcarni et al. 1998). The fire- 
ball observations immediately after the burst, when the tem- 
perature and density are at maximum, are of great interest. Un- 



fortunately, it has been impossible to observe the afterglows 
in the radio, optical, or X-rays earlier than approximately 10 
hours after the burst. 

Some earlier observations indicated that the afterglows 
could immediately follow some GRB. There were detections 
of X-ray emission lasting for tens of seconds after the main 
burst was finished in gamma-rays (Sunyaev et al. 1990, Mu- 
rakami et al. 1991, Terekhov et al. 1993, Sazonov et al. 1998). 
PVO observatory observed a faint gamma-ray emission over 
~ 1000 s after flie long, 200 s, burst GRB 840304 (Klebesadel 
1992). The presence of slowly fading soft gamma-ray (100- 
500 keV) emission was found in about ^ 10% of bursts de- 
tected by GRANAT/PHEBUS (Tkachenko et al. 1995). Hard 
gamma-ray photons (0.2-10 GeV) were detected during 1.5 
hours after GRB 940217 by EGRET telescope (Hurley et al. 
1994). 

We present here a detailed analysis of the GRB 920723 
light curve, which reveals a soft gamma-ray afterglow with flux 



decaying as a power law 
main burst. 



2. Observations 
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during at least 1000 s after the 



SIGMA is the coded-mask telescope with a 15' angular reso- 
lution operating in the 35-1300 keV energy band (Paul et al. 
1991). Typicafly, SIGMA performs uninteiTupted 20-30h ob- 
servations, during which the telescope pointing is maintained 
with a 30' accuracy (but known to within 15"). Although the 
telescope field of view is only 11.4° x 10.5° (FWHM), some 
fraction of gamma-rays from sources closer than 35° to the 
pointing direction reaches the detector through the gaps in 
the passive shield and produces arc-shaped images. This "sec- 
ondary optics" (or "sidelobes") was described in detail by 
Claret et al. (1994a) along with the appropriate analysis tech- 
niques. 
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Fig. 1. GRB 920723 light curve with time resolution ^ 0.1 in 
35-300 keV band. The reference time is at the burst trigger. 
Vertical dotted lines represent the reference times of the light 
curve in logarithmic coordinates (Fig. 2 and 4). 



GRB 920723 was observed by SIGMA through the sec- 
ondary optics. An ^1° localization was obtained from this 
observation (Claret et al. 1994b). GRB 920723 is one of 
the brightest bursts observed by GRANAT instruments and 
the brightest detected by SIGMA. The burst was triggered at 
20''03'"08\3 UT and lasted for about 6 s. The WATCH all- 
sky monitor provided a 0.2° localization (Sazonov 1998) and 
observed the fading X-ray emission in the 8-20 keV band dur- 
ing more than 40 s after the main burst (Terekhov et al. 1993). 
PHEBUS measured the peak burst flux 5 x 10"^ erg s~' cm~^ 
and fluence 1.4 x 10"^ erg cm"^ in the 100-500 keV energy 
band (Terekhov et al. 1995). 

The SIGMA data allows the measurements of the burst 
light curve with better than 0.1 s time resolution (depending 
on flux) during 7.5 s after the trigger. In addition, the count 
rate in four wide energy bands (35-70, 70-150, 150-300, 300- 
600 keV) is recorded with the 4 s time resolution over the entire 
observation. With these data, it is possibile to study the burst 
emission long after the trigger. Below, we use only the first 
three energy channels, because the last one is plagued by low 
sensitivity. The peak burst count rate in the 35-300 keV band 
was 7900 cnt s"', much higher than the average background 
rate 310 cnt s"^ 

During the observation of July 23, 1992, SIGMA was 
pointed to Her X-1. The pulsar was in eclipse and was not 
detected. The 3cr upper limit on its 35-70 keV flux averaged 
over the observation, was 0.25 cnt s"'. The Her X-1 spectrum 
is known to be very soft (the 20-100 keV photon index is -4.4), 
and therefore its flux is negligible above 70 keV. The pulsar was 
in eclipse between 12000 s before the burst and 9000 s after 
the burst. Therefore, it could not cause any significant variabil- 
ity of the SIGMA count rate during the reported observation. 
No other known bright sources were visible through either pri- 
mary or secondary optics of the SIGMA telescope. GRANAT 
operates on the high apogee orbit and, during the observation, 
was not influenced by the Earth radiation belts or other mag- 
netospheric anomalies (such as the South Atlantic Anomaly). 




Fig. 2. The background-subtracted burst light curve. Zero time 
is at the burst trigger. The filled squares show the count rate 
measured with the 4 s resolution, and the open squares — the 
data with good time resolution. The horizontal dotted line rep- 
resents a 95% upper limit on the possible internal background 
variations on the 300 s time scale (note that Poisson variations 
are already included in the error bars). 



As a result, the SIGMA background usually does not show any 
significant variations on the time scales shorter than ^ 10^ s. 
Therefore, it can be accurately modeled by a low degree poly- 
nomial. 



3. Results 

Usually, sources contribute only a small fraction of the to- 
tal SIGMA count rate. Therefore, the correct background sub- 
traction is vital for the source variability studies. We modeled 
the background using Chebyshev polynomials. A complete de- 
scription of the SIGMA background subtraction techniques is 
presented elsewhere (Burenin et al. 1999); this analysis has 
shown that the background variations around the subtracted 
value in excess of Poisson noise are smaller than 0.6 cnt s"' 
on the 300 s time scale (on the 95% significance level). 

Figure 1 shows the burst light curve in the 35-300 keV 
band. There is a small peak in burst light curve just after the 
trigger. Within the first second after the trigger, the burst flux 
rose rapidly. Over the next five seconds, it remained at approx- 
imately the same level, showing a strong variability at all re- 
solved time scales. At approximately 6 s after the trigger, the 
flux started to decline rapidly. 

Figure 2 shows the burst light curve in logarithmic coor- 
dinates of both time and flux. The shape of the light curve in 
these coordinates strongly depends on the choice of the refer- 
ence time. In Fig 2, the reference time is chosen at the moment 
of the burst trigger. There appears to be a power-law decay of 



Fig. 3. Time history of the effective spectral index in the 8- 
200 keV energy band. Zero time is at the burst trigger. Vertical 
dotted line represents the moment f = 6 s from the trigger, when 
the burst flux began the gradual decline. 
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Fig. 4. Same as Fig. 2, but the reference time was set at 6 s 
after the trigger. The main burst is not shown here because it is 
at f < with this choice of reference time. 

flux starting at 10-20 s after the trigger. This behavior is con- 
sistent with GRB entering the stage of self-similar fireball ex- 
pansion soon after the main burst. It should be emphasized that 
the self-similar behavior is expected only on time scales much 
larger than those of the main energy release. Therefore, we used 
the data in the 20-1000 s time interval for the power-law mod- 
eling of the light curve. 

The solid line in Fig. 2 shows the power law fit in the time 
interval 20-1000 s; the dash-dotted line shows the exponen- 
tial fit in the same interval. The reduced x^ is 1-5 and 6.5 (4 
dof) for the power law and exponential models, respectively. 
The power law adequately describes the data and results in a 
better fit than the exponent. The best fit power law index is 



-0.69 ±0.17 (Ax^ = 2.7 for the index -1). This power law tail 
contains at least r^ 20% of the main burst fluence. Applying the 
same procedure to the three wide energy channels separately, 
we obtained the power law indices -1 .29 ± 0.55, -0.64 ±0.19, 
and -0.41 ± 0.6 in the 35-70, 70-150, and 150-300 keV en- 
ergy bands, respectively. Interestingly, the extrapolation of the 
power law (shown by the dotted line in the Fig. 3) points to the 
small peak near the beginning of the main burst (Fig. 1). 

The spectral evolution of the burst flux can be character- 
ized by the ratio of the 8-20 keV flux measured by WATCH 
(Terekhov et al. 1993, Sazonov et al. 1998) and the SIGMA 
flux in the 75-200 keV band. Standard SIGMA calibration does 
not apply to the secondary optics flux, so we used the PHE- 
BUS measurement of the main burst flux to find the conversion 
coefficient between SIGMA counts and flux. Figure 3 shows 
the observed ratio of 8-20 and 75-200 keV fluxes expressed 
in terms of equivalent spectral index in the 8-200 keV energy 
range (i.e., the spectral index of an /v ^ z^"" spectrum with the 
same hardness ratio as observed). The time behavior of a in- 
dicates that the afterglow spectrum is significantly softer than 
that of the main burst. During the main burst (in the 0-6 s time 
interval), the usual hard-to-soft evolution of the GRB spectrum 
is observed (e.g. Ford et al. 1995). Near the start of the grad- 
ual decline of the burst flux, 6 s after the trigger (shown by the 
vertical dotted line in Fig. 1 and 3), a changes abruptly from 
«0.3to«l. 

It is interesting to examine the light curve with the reference 
time chosen at f = 6 s after the trigger because this moment can 
be singled out in both flux and spectral history of the burst; 
the result is presented in Fig. 4. With this choice of zero time, 
the data in the 0.01-20 s time interval lies on the extrapolation 
of the power law fit from the 20-1000 s time interval. Adding 
these data to the fit results in the power law index -0.70 ± 0.03. 
Note that in this case, the power law flux decay is observed over 
approximately four orders of magnitude of time. 

4. Discussion 

We presented a high sensitivity observation of the GRB 920723 
light curve in the soft gamma-ray band. The stable background 
of SIGMA allows detection of the burst emission on the level 
of better than 1/1000 of the peak intensity. A similar analysis 
would be complicated with BATSE because the background is 
less stable and because the source is eclipsed by Earth every 
several thousands seconds. We were able to detect the burst af- 
terglow extending up to ~ 1000 s after the main burst. There 
is a continuous transition of the main burst to its power law af- 
terglow (Fig. 1 and 2). The afterglow spectrum is significantly 
softer than that of the main burst. An abrupt change in the burst 
spectrum occurs at approximately the same moment when the 
power law decay of flux seems to start, at f « 6 s after the trig- 
ger 

The behavior of GRB afterglows in the lower energy bands 
and at f > 3 X 10'* s can be explained by the synchrotron emis- 
sion of electrons accelerated in external shocks generated by 
relativistically expanding fireball colliding with the interstel- 
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lar medium (e.g. Wijers et al. 1997, Waxman 1997, Sari et al. 
1998, Wijers and Galama 1998). In the framework of this 
model, the spectral flux at the observed frequency v is given 
by Fi, oc v~°T^^ , where a and ji are constant and depend only 
on the spectral index of electrons on sufficiently late stages of 
flie fireball evolution. For GRB 920723 we obtain a = 1 ± 0.2 
and /9 = 0.69 ±0.17. Both the spectrum and the light curve of 
GRB 920723 seem to be considerably flatter than that of X-ray 
afterglows observed for other gamma-ray bursts at f > 3 x lO'* 
s — a= 1.4-1.7 and/3= 1.1-1.6 (e.g., in't Zand et al. 1998, 
Piro et al. 1998, Nicastro et al. 1998). Furthermore, the flux 
decay observed by SIGMA is flatter than r' (at ^ 90% con- 
fidence), i.e. the total flux diverges if extrapolated to f ^ oo. 
This suggests that the afterglow light curve should steepen at 
some moment during or after the SIGMA observation. 

In the relativistic fireball model, the afterglow light curve 
and the energy spectrum should steepen simultaneously at the 
moment f„, when the maximum in the electron spectrum, E„„ 
passes through the SIGMA bandpass (we assume below that 
t,n corresponds to E„ = 100 keV). At later stages of the fireball 
evolution, indices a and [3 do not change with time. Since the 
light curve steepening after 100-lOOOs is required by our data 
at ~ 90% confidence, it can be suggested that f,„ > 100- 1000 s. 
Also, if indeed t,„ is > 100 s, our flat spectrum may become 
consistent with the parameters of X-ray afterglows (see above) 
because the spectrum softens after t,„. In the adiabatic fire- 
ball, t„t can be estimated as f„, « 140 e^ Se £53 s, where 
£53 X 10^^^ ergs is the total energy release, and Se <l and Eb <l 
are the fractions of the electron and magnetic field energy in the 
total shock energy, respectively (Sari et al. 1998). The value of 
t,„ not much less than 100 s would not be strongly inconsis- 
tent with the SIGMA data. However, using Se and Eb estimates 
from the parameters of radio, optical, and X-ray afterglows of 
GRB 970508 and GRB 971214 at t - lO^-lOS (Wijers & 
Galama 1998), we obtain t„, ^ 3 s, which does seem to contra- 
dict our data at ^ 90% confidence level. This may indicate a 
large diversity of the fireball parameters in different bursts or 
some problems of a simple model of a spherically symmetric 
fireball in explaining the early stages of the gamma-ray burst 
afterglows. 

SIGMA data provides the first convincing observation of 
the power law afterglow in the soft gamma-rays and immedi- 
ately after the burst. A very important issue is whether such af- 
terglows are common. A preliminary analysis of other SIGMA 
bursts revealed no other convincing afterglows, primarily be- 
cause of the faintness of other bursts; on the basis of SIGMA 
data alone, we cannot rule out that the soft gamma-ray after- 
glow is a common phenomenon. A preliminary analysis of 
the PHEBUS data confirms the detection of the afterglow in 
GRB 920723 and reveals a similar afterglow in GRB 910402 
(Tkachenko et al. 1998). The results of our systematic search 
for soft gamma-ray afterglows in the GRANAT data will be 
presented in the future. 
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Appendix A: Background modeling 

(This section does not appear in the Journal version) 

We start with fitting the background by the Chebyshev 
polynomial. We exclude the data between 1000 s before the 
burst and 4000 s after the burst. We increase the order of the 
polynomial until the F-test indicates that no additional pow- 
ers of t are necessary; we find that the polynomials of the sec- 
ond and third orders are required to describe the background in 
spectral channels 1&2 and 3, respectively. 

Since the fit is made for the entire observation, the uncer- 
tainty of the fit value, arising from statistical uncertainties in 
the polynomial coefficients, is negligible in any small part of 
the observation — at least compared to the statistical error of 



flux in that part. So, we do not consider the fit uncertainties any 
further. 

The deviations of the data from the fit should ideaUy be 
Poissonian. However, we cannot exclude apriori the existence 
of internal background variations on any time scale. So, we 
want to place an upper limit on such internal variations. We 
proceed as follows. 

a) We choose the time scale 300 s for this study because 
this is the width of time bins near 1000 s in Fig. 2 and 3. 

b) We average the observed flux in the 300 s bins and make 
a histogram of deviations from the fit expressed in units of Pois- 
son error in this bin. In the absence of internal variations, this 
histogram should be consistent with the Gaussian with zero 
mean and standard deviation = 1 (i.e., with the normal distri- 
bution). 

c) We do find that the histogram can be described by normal 
distribution. This shows that there are no biases in the back- 
ground determination and that the internal background varia- 
tions on the 300s time scale are small. To set the upper limit, 
we fit the width of the distribution and then convert the upper 
limit of the width into the corresponding count rate, assuming 
that internal variations (if any) and Poisson noise were added 
in quadrature. 

This technique results in a 95% limit of 0.6 cnt/s for internal 
background variations on the 300 s time scale. 



